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Abstract. We undertake a non-perturbative study of the evolution of the 
“gravitational entropy” proposed by Clifton, Ellis and Tavakol (GET) on local 
expanding cosmic CDM voids of ^ 50—100 Mpe size described as spherical under-dense 
regions with negative spatial curvature, whose dynamics is determined by Lemaitre- 
Tolman-Bondi (LTB) dust models asymptotic to three different types of FLRW 
background: ACDM, Einstein de Sitter and “open” FLRW with A = 0 and negative 
spatial curvature. By assuming generic nearly spatially flat and linear initial conditions 
at the last scattering time, we examine analytically and numerically the GET entropy 
evolution into a fully non-linear regime in our present cosmic time and beyond. Both 
analytic and numerical analysis reveal that the late time GET entropy growth is 
determined by the amplitude of initial fluctuations of spatial curvature at the last 
scattering time. This entropy growth decays to zero in the late asymptotic time range 
for all voids, but at a faster rate in voids with ACDM and open FLRW backgrounds. 
However, only for voids in a ACDM background this suppression is sufficiently rapid 
for the GET entropy itself to reach a terminal equilibrium (or “saturation”) value. 
The GET gravitational temperature vanishes asymptotically if A = 0 and becomes 
asymptotically proportional to A for voids in a ACDM background. In the linear 
regime of the LTB evolution our results coincide, qualitatively and quantitatively, 
with previous results based on linear perturbation theory. 
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One of the long standing open issues in current theoretical physics is the proper dehnition 
of a “gravitational” entropy, Sg^, that determines the directionality of the gravitational 
interaction at all scales and energies. Hence, it is a different (though possibly related) 
notion from that of the thermal entropy of the held sources (hydrodynamical or non- 
collisional) or the holographic black hole entropies. Research on this issue has produced 
various self-consistent proposals, from Penrose’s old idea of the “arrow of time” [1] 
(see also 0 ). to the proposals by Clifton, Ellis and Tavakol [3] (GET) and Hosoya 
and Buchert (HB) |3] (see comprehensive study of these two proposals in [S] and other 
related research in 0 [3 El El [la E]). 

In a cosmological context the FLRW models emerge in the various proposals of this 
entropy as a global “equilibrium state” (zero entropy growth Sg^ = 0 for all fundamental 
observers), thus characterising the deviation from homogeneity and isotropy (inherent in 
structure formation and cosmic expansion) as irreversible processes whose viability can 
be verihed by means of a theoretically consistent test: demanding that astrophysical and 
cosmological solutions of Einstein’s equations (whether exact, numerical, perturbative 
or post-Newtonian) fulhll the conditions for gravitational entropy growth Sg^ > 0. 

Considering the empirical success of the ACDM paradigm [Ellis], as well as the 
fact that numerical simulations m reveal a CDM density distribution made of a web of 
local cosmic voids of ~ 50 — 100 Mpc size, we believe that these void structures within 
a concordance cosmology background provide an ideal non-perturbative framework to 
probe the evolution of the GET entropy proposal, thus complementing a recent study by 
Marozzi et al HDl and earlier work mm on perturbative studies appropriate for large 
scale structure. A convenient non-perturbative approach follows by considering idealised 
spherically symmetric CDM voids described by Lemaitre-Tolman-Bondi (LTB) dust 
models with negative spatial curvature [151 HB]- In particular, we aim at looking at 
the effect of a nonzero cosmological constant by comparing the GET entropy growth 
predictions for voids whose asymptotic FLRW background is a ACDM model (hence 
A > 0) and an FLRW background model with A = 0 {i.e. spatially flat Einstein de Sitter 
or open FLRW with negative spatial curvature). While voids in an FLRW background 
with A = 0 can be ruled out on the grounds of compliance with observations, it is 
nevertheless important from a theoretical point of view to examine the suppression 
effect of negative spatial curvature on the growth of the gravitational entropy. 

We have shown in previous work 00 that, as long as the decaying mode is 
subdominant, a growing GET entropy Sg,. > 0 is a generic feature of LTB models, 
whether ever expanding or collapsing and irrespective of the void or over-density prohles 
or the existence of an FLRW background. We also showed that for ever expanding 
models the rate of growth Sg^ is eventually suppressed for all fundamental observers, 
with a decay Sg^ —)■ 0 in the asymptotic time range. We extend this work by computing, 
through analytic asymptotic expressions, the specihc functional form of the entropy 
growth suppression. The result is a rapid exponential decay of Sg^ for voids in a ACDM 
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background and a slow logarithmic decay for voids in FLRW backgrounds with A = 0 
(Einstein de Sitter or open FLRW). As a consequence, the GET entropy itself, Sg,, which 
follows from a suitable integral of the time rate Sg^ along the worldlines of fundamental 
observers, only reaches a finite (position dependent) asymptotic terminal “saturation” 
value for voids in a ACDM background, diverging logarithmically for voids in FLRW 
backgrounds with A = 0. 

In order to complement the analytic work described above and to obtain 
quantitative results, we examine (numerically) the evolution of the GET entropy Sg^ 
and its rate of growth Sg^ for 50-100 Mpc size voids with the three FLRW backgrounds 
(AGDM, open FLRW and Einstein de Sitter), all evolving from the same plausible 
generic linear initial conditions (at the last scattering surface) and with a background 
value for Hq compatible with cosmic age constraints. We obtain results consistent with 
the asymptotic analytic expressions. We find that the rate of growth Sg^ is determined 
by the amplitude of initial fluctuations of the negative spatial curvature and is basically 
insensitive to initial fluctuations of the GDM density. The growth suppression proceeds 
very slowly in the void in an EdS background, so Sg^ is still growing at a fast rate in our 
cosmic time t^. For the voids in the AGDM and open FLRW backgrounds a noticeable 
suppression of Sg^ is already present at relatively early times t ~ to/2, proceeding at a 
faster rate for the AGDM background, though s^, itself is still larger att = to than in the 
void with an open FLRW background. However, for cosmic times beyond to the growth 
suppression becomes substantially larger in the void with a AGDM background, with Sg,. 
reaching a terminal value at ~ 2to, while Sg^ keeps growing at a very slow logarithmic 
rate in the void with an open FLRW background. Hence, negative spatial curvature 
does suppress the GET entropy growth, but not at a sufficiently fast rate to allow for 
the asymptotic convergence of Sg^. 

The analytic and numeric results summarised above are qualitatively analogous 
to and consistent with those obtained in the linear perturbative analysis of [10], who 
find that the GET entropy growth is suppressed (he. “saturated”) only for dust 
perturbations in a AGDM background, with this entropy diverging asymptotically for 
their “GDM Universe” (an Einstein de Sitter background), as these authors did not 
examine the effects of negative spatial curvature that would arise from an open FLRW 
background with A = 0. Since the exact fluctuations that we used to examine the 
non-perturbative evolution of s^^ and Sg,. reduce to standard cosmological perturbations 
in the linear regime mi, we examine the entropy growth in a linear regime that is 
fully consistent with spherical dust perturbations in the isochronous comoving gauge. 
As expected, the rate of growth in this regime is much smaller than in the non-linear 
regime. 

The GET “gravitational temperature” associated with the terminal equilibrium 
value of the GET entropy is in all cases proportional to the Hubble expansion scalar, 
hence we can identify for local voids in a AGDM background an asymptotic equilibrium 
temperature proportional to \/A, whereas for the open FLRW background this terminal 
temperature is zero. 
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The paper is organised as follows. In section we briefly present LTB models 
described in terms of covariant qnasi-local scalars (the “q-scalars”). We introdnce 
in section the GET entropy and its associated gravitational temperatnre. Initial 
conditions at the last scattering time and dimensionless variables are introdnced in 
section An analytic stndy of the rate of growth of the GET entropy and temperatnre 
in the asymptotic time regime is provided in section while in section we examine 
the snppression of the GET entropy integral and its asymptotic time limits. In section 
[^we examine throngh nnmerical examples of void strnctnres the evolntion of the GET 
entropy and its rate of growth, with a detailed comparison with a linear pertnrbative 
approach in section]^ Finally, we snmmarise in section [^onr resnlts and conclusions. 


2. LTB dust models. 


A convenient parametrization for the LTB dust models is furnished by the following 
FLRW-like metric (we use units with G = c = 1)\ 

" r2 


ds" = 


-dt^ + 


I-1C. 


QLS' 


j dr^ + (d-d^ + sin^ 'dd(p^) 


r a 


T = r(f,r) = 1 H -, a = — 


a = a(t,r) satishes 

3 V a / 


da 
dr ’ 


ICn 


( 1 ) 

( 2 ) 

( 3 ) 


da 

qhS, , 

where the r-dependent functions /C^ls = lCg{Gs,r) and p^Ls = PqitLs,^) are dehned 
further ahead in (Q. The subindex ls will denote henceforth evaluation at the last 
scattering time f = ^ls ~ 4 x 10® ys. The radial coordinate has been chosen so that 
Ols = Tls = 1 hold for all r. 

The standard procedure in most applications of LTB models is to use the solutions 
(analytic or numerical) of the Friedman equation (|^ (in various parametrisations) to 
determine a and T and then compute all relevant quantities (usually the present cosmic 
time t = to is used as reference hducial time). An alternative approach is to consider 
the covariant fluid flow scalars of the models: the energy density p, the Hubble scalar 
T-L = 0/3 with 0 = h’fWhU°‘ and the spatial curvature K, = ^^^TZ/Q with the Ricci 
scalar of hypersurfaces of constant t (orthogonal to m“). However, the “quasi-local” 
scalars (to be denoted by “q-scalars”) and their fluctuations provide a set of equivalent 
and useful covariant scalars that are specially suited to probe the GET entropy [3] 
and reduce in a linear regime to gauge invariant spherical dust perturbations in the 
comoving isochronous gauge Ha- These scalars (which will be denoted by the subindex 
q) are dehned as [f] 

1C. 


r: AR^R'df 

A — do_ 

/o R^'dr 


Pq = 


PqhS 


'Hq = -, 
a 


lCq = 


■.'IJLS 


( 4 ) 


f The q-scalars Aq are closely related to weighted proper volume averages of the scalars A. They have 
been very useful to study theoretical properties of LTB models [Miiiiiiiiniiii] and Szekeres models 
[22| . Their properties are comprehensively examined in these references. 
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where R = ar is the area distance, and they are related to the “standard” covariant 
scalars A = p, R, )C throngh their flnctnations 

r 1 

= A-Ag = ^ = —J A' R^fdf. (5) 

The q-scalars and their flnctnations satisfy the following algebraic constraints 

+ 2W,D‘«> = - dot. (6) 

The shear and electric Weyl tensors are nicely expressed in terms of the flnctnations of 
the density and Hnbble scalar: 




K = = T 

a 0 a 




(7) 


where ^2 is the conformal Newman-Penrose invariant and — Suan^ = 

diag[0, —2,1,1] (with a nnit vector orthogonal to and to the orbits of SO(3)). 
As shown in previons literatnre [3 na [la Eoi EB Ea Ea], the q-scalars and their 
exact flnctnations provide a representation of covariant scalars that fnlly determines 
the dynamics of LTB models (either analyticaly or nnmerically). 


3. The GET gravitational entropy in LTB models. 


The GET gravitational entropy is defined from an “effective” energy momentnm 
tensor for a snitable geometric held associated with the Weyl tensor (the “free” 
gravitational held [§]). For LTB models (which are Petrov type D spacetimes or 
“Conlomb-like” fields), this tensor takes the form [3]: 

^ + n“?. (8) 

where is the 4-velocity of the matter sonrce and the “gravitational” density, pressnre, 
viscosity and heat flnx are given in terms of the conformal invariant \k 2 

(directly related to the density flnctnation by ([^) as 


Q/jT 

STTPg, = —a 


dvr 


87rn“(’ = — 
o 

where = x°‘x'^ + y^y’^ 


ab 


Pgr = = 0, (9) 

+ u°‘u^ for the canonical orthonormal tetrad 
{m“, n“, ?/“} associated with ([^ and a is a constant to get the right nnits. By 

setting np a formal analogy between in ([^ and the energy-momentnm tensor of 
a “reference” dissipative hydrodynamical sonrce with 4-velocity [25], GET obtain a 
Gibbs eqnation analogne for the gravitational entropy growth, which for LTB models 
takes the from (see details in [3] and 0 ): 


(pgr^) ^ ^ab 


gr ' 


P 




E 


ab 


a 


27ra d 


^(|D<qa3r) 


= - 27rap, a^PD^^) 


D(p) 


( 10 ) 


§ GET obtain the second order “effective” energy-momentum tensor through an irreducible 
algebraic decomposition of the Bell-Robinson tensor, the only fully symmetric divergence-free tensor 
that can be constructed from the Weyl tensor. See comprehensive discussion in [3|. 
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where we used ([i]), (36a)-(36d), (36e), the tensors aah^ follow from Q, Id = G" = a^T 
is the local volume defined hj £/i = T-L and the “gravitational” temperature was given 
by GET in [3] as 


T = 


\n + aabn'^n^\ _ \nq + 3I)^^^\ 


( 11 ) 


2n 271 

on the grounds that it reduces to the semi-classical Unruh and Hawking temperatures in 
the appropriate limits (see discussion in [3]). The terms inside the brackets in the right 
hand side of (10) describe the “effective” dissipation produced by the GET entropy 


through the geometric fluid associated with ([^ (notice that there is no exchange of 
energy or momentum between this formal fluid and the LTB dust source). Since we 
have aab = 0 for FLRW models, the latter define a global “gravitational” equilibrium 
state characterized by = \'H\/{27i) > 0 and Sg^ = 0 for all t and all fundamental 
observers. 

Assuming that T > 0 holds to avoid shell crossing singularities [iniEi], we obtain 


from (10) the necessary and sufficient condition for entropy growth along the worldlines 


of fundamental observers [5]: 

s„ > 0 » < 0. 


( 12 ) 


which takes the form of a negative correlation between fluctuations of the energy density 
and the Hubble scalar with respect to their q-scalars. 


As argued in |5], the GET entropy growth condition (10) is only the time component 
(he. component projected to the 4-velocity) of the one-form 


Tg, dsg,. = -7iad{p^^V). 
o 


(13) 


whose integrability conditions for LTB models place the following restriction on the 
radial dependence of s^, (see details in |5]): 

{Tiy . • h ; + 3 ( d (^))' . 


s = 


^2y 
' grl 

iTl) 


Sa-r - 


Ha + 3D 


in) 


(14) 


which must be considered when looking at the integration of dSg^ in (13) (see section 

§. 


4. Initial conditions and dimensionless variables. 


Since we will consider initial conditions at the last scattering time t = tj^s, it is useful to 
normalise all initial value functions with respect to the Hubble length scale of the 
FLRW background at this time (an over-bar will denote henceforth FLRW background 
quantities). In particular, the initial GDM density and spatial curvature (standard 
scalars not q-scalars) and cosmological constant can be given by 

_ _ ^(fc) 


2fiLs = TT^ = + rniC^^Hr] 




\ — o(A _ o(A) ^0 

^ ~ ^^0 TTO 


= 1.37 X 10 


-9 


^LS — 


(Hq^^ = 0.73, Hqc = 69km/s Mpc), 


772 


= + hy^\r), 


(15 a) 


(156) 
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where 77 ls, and — 1 characterise the FLRW background at 

Gs, the smooth functions e^'^\r), e^^)(r) are selected to fulfil regularity conditions (see 
[20]) and to satisfy = e^^^(O) = 1 and —)■ 0 as r —)■ oo, while the constants 

mi, ki comply with (|mi|, |fci|) ^ 1. The q-scalars and fluctuations associated to /Xls 
and Kls are 


/Tijls — 


47rp, 


gLS 


Q 


(m) 

LS 


ms 


+ mic, 


(m) 


l^qLS — 


8'^^gLS _ q(^) 

ms "" 


+ hei’‘\ 


h - 

ff'qhS — 


h 

•^qhS ~ 


H,s 

1 + mic 


— [2p-qLs ~ f^qLs + ~ 4" 2mie^™'^ — 

M _ b.Ak) 


1/2 


2^^ 


— 

^^qhS — 


h‘2 

^qGS 


n 


(A) 

LS ? 


= "ms 




-G = -li 

gr-gija; “ g'^lJLS) -^LS 2 "'(JLS 

where we applied (0)-@ to ( |l5a| -( [l^ and are obtained from from 

the integral (|^ evaluated at Rs (he. with R = r, since Ols = Tls = 1). 

For further studying the evolution of the GET entropy for the initial conditions 
given above, we introduce the dimensionless time and radial coordinate normalised with 
respect to the Hubble scale 


d["s^ = m. = 


21^ 

2h, 




(16 a) 
(166) 
(16c) 
(16d) 


gLS 


T = Hr.Ct - 1 


3 \ t 


X = H^sr = 


2r 

3ctr, 


(17) 


SO that the initial time slice t = Rg is given by r = 0 and present day cosmic time 
to ~ 13.7 Gys is Tq ~ 2.2 x lO'^, while y = 1 corresponds to the Hubble radius 
r = 1 /Hls = 3ctLs/2 ~ 0.19 Mpc. 

Regarding the entropy growth, we recall that the term 27ra in (10) has units 
of entropy times length (after writing 27rpg as 2T\Gpql(?'). Hence, considering that 
Pqof = PgLs, the components of the GET entropy one-form (13) in (10) and (14) can be 
written in terms of dimensionless quantities as 

9S„ 3 row DO-' 

dT ~ ( 18 ) 


|hg + 3D('*)| |D(^)|’ 


Qx 


h-lJLS 


hq + 


-I ,x 


X 






(19) 


where Sg,. = Sg,./cT, with a being a suitable quantity with entropy units (its theoretical 
role is discussed in section and 


Yin — 


77t. 


= h-hn = 


d(«) 

m' 


=P-Pq = 


47rD(^) 

3Hh ' 


( 20 ) 


generalise for all r the initial value functions hgLs, Dj/g and D/^g . The two components of 
dSgr in (18) and (19) must be considered separately when looking at along comoving 
observers {ds^Jdr for y constant) and along radial rays of hypersurfaces orthogonal to 
the 4-velocity {ds^Jdx for r constant). 
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5. Asymptotic entropy growth and temperatnre. 


We proved in [HI that ds^Jdr > 0 holds in the time asymptotic regime of fundamental 
observers in generic LTB models with negative spatial curvature in the case A = 0. We 
review this result here and prove that it also holds for the case A > 0. We also examine 
the asymptotic evolution of the gravitational temperature Tg,. For this purpose we need 
the time asymptotic forms of the metric functions a and F, which can be obtained from 
the asymptotic range a 3> 1 (or f/^Ls 3> 1) of the quadrature of the Friedman equation 


(§ 


Vt 




(m) 

gLS 




(A) 


nl/2^ 


where the Omega factors are 


— 


2/^ijLS 

""gLS 


— 

""" - hi 


A 


gLS 


o(0 _ _ Q(m) 

— 7 2 — '^^g'LS 

^qLS 


+ Q 


(A) 

gLS 


1 , 


( 21 ) 


( 22 ) 


and the lower bound of the integral corresponds to t = ^ls (or r = 0, a = Ols = 1). For 
ever expanding dust layers (see conditions for this when A > 0 in |23]) the late time 
forms for the metric functions a and F are 


exp 


qGS 


1 I QaC^) “T iigrLSj^LS ^LS 

1 + oOls 0 




1 _ h r 
1 iigLS ' ? 


1 + ^4? + 


30, 


(m) 

gLS 


Av) 

Ogs 


SAG 

2"ls 


In a 


1 - 0 , 


(m) 

gLS 


(A > 0), 

(A = 0), 


where we have introduced the relative initial fluctuations 


AG _ 

Gjo — 


hgLS 


AG _ 




= 




gLS 


hqLS 


Oim.) 

“gLS r(At) 


qIG 

^^gLS ^(k) 


(23) 


(24) 


Inserting the forms (23) into (|^ and (|^ (expressed in terms of dimensionless quantities) 
and considering only leading terms in 1/a leads after a long algebraic manipulation to 
the time asymptotic form of the product of fluctuations 


qIA j3(a) 


'^gLS^^gLS 


q(^)AG _ o(ny"o I n 

i^gLS^LS ^l^i^OLS ^ 


){m) 

^gLS 


{A)AeW 

gLSj^LS 


_ oAG 

^LS ’J^LS 


4a/0^^^ 

(^LS 


gLS 

2h 


(l + 35S))2a5 


LS 


')(5^^^-3<5S)) 




(A > 0), 


(25 a) 


D{G _ 


[(2 + AG& - 63S> 


In a 


4i/l-!jg 




o(™) 

^ ^gLS 

35[s^ 


1 + ) - 1 - 35, 


In a 


(h) 

LS 


1 - O^”"^ 

A ^ ^gLS 


a(/^) _ 3 r(^) 
^LS 2 ^ls 


(A = 0), 


(256) 
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where we used the approximations k, 1, k 


gLS 




•C 1. The expansions 


(25a)-(256) determine the late time sign of ds^Jdr through condition (12) for arbitrary 
radial profiles of initial value functions specified at ^ls- 

Considering the fact that —)■ 0 holds as r —)■ cxd [a ED, the gravitational 

temperature associated to the asymptotic equilibrium state defined by ds^Jdr —)■ 0 is 
given by —)■ hi„f/27r where is the limit of h = as r —)■ oo. We look at the 

cases A = 0 and A > 0 separately below. 


5.0.1. The case A = 0. As proven in |5], equations (12) and (256) confirm that the 
fulfilment of ds^Jdr > 0 in the asymptotic time range of fundamental observers is a 
generic feature in hyperbolic models with negative spatial curvature (0 < < 1 

implies Kg < 0 for all f), irrespective of the nature of the density profile at f = ^ls or of 
the existence of a decaying mode which is already subdominant at ^ls and is completely 
negligible for late times mm- We have Anf = 0, hence Tg^ —)■ 0, and thus the asymptotic 
equilibrium state is characterised by zero gravitational temperature. 


5.0.2. The case A > 0. Considering a ACDM background with = 0, together with 
the relation between the signs of the radial gradients and the fluctuations in (|^, and 
assuming monotonic (but otherwise unrestricted) radial density and spatial curvature 
profiles at the initial slice we have: 

• Initial void profile: ^ 0) hence > 0 (density minimum at the void centre), 

but AgLs < 0 and thus < 0 (maximal positive Hubble expansion at the centre 
of the void). 

• Initial over-density profile: /ig^g < 0, hence < 0 (density maximum at the 
clump centre), but hg^g > 0, hence (f^g^ > 0 (minimal positive Hubble expansion at 
the clump centre). 


and thus it is straightforward to show from (12) and (256) that ds^Jdr > 0 holds in 
the asymptotic time range. While the numerical examples in |23] show this result, we 
prove it here analytically. Since we have 


27rT'gr —)■ 


= \/yA. 


(26) 


the CET entropy provides (irrespective of initial conditions) a nice interpretation of A 
as proportional to the asymptotic gravitational temperature of the terminal equilibrium 
state. 

6. The suppression and “saturation” of the CET entropy. 


The correlations in (25a) and (256) suggest a dramatic asymptotic suppression of the 


CET entropy growth for fundamental observers when A > 0. We still need to verify this 
suppression (and the asymptotic suppression in the radial direction) from the general 
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expressions (18) and (19). This information will allow ns to infer the general asymptotic 
behavionr (for fnndamental observers) of itself. 


6.1. Entropy suppression for fundamental observers. 


Considering the fnll form of ds^^Jdr from (18), nsing the expansions (23) and the 


prodncts of the asymptotic form of the rate of decay of entropy growth for 

fnndamental observers (y constant) is given by: 


oA) 

ds., 3^ Vs 

° rs_/ 

f _i_ 

ar 8 



ds., 3^ Vs 

° rs-/ 

qQ{^)rA{^) 

^qhS ■V'lS ^^LS-V'lS 

In a 

r«s_; 

dr 4 


a 




.-2' 


8U 


(A)' 

qhS 




(A > 0), 


9D|,"s^ 


Inr 


qhS) 


4(a 


G) \2 

qLS) 


T 


(A = 0), 


(27) 


(28) 


which reveals a decay to zero of entropy growth rate for fnndamental observers located 
in voids in all FLRW backgronnds (irrespective of whether A > 0 or A = 0). However, 
this decay occurs at a fast exponential rate in voids with a ACDM background, while 
it is a slow logarithmic decay Int/t for voids in any FLRW background with A = 0. In 
order to explore the consequences of these decay rates, we introduce the entropy integral 
for fundamental observers {i.e. the line integral of the one-form ds projected along the 
integral curves of the 4-velocity) 


S^AT 


{t)\x = / (dsgr ■ u)dr = f){x) + 


r 

ds^,. 

lo 

. dr _ 


dr. 


(29) 


^=const 


where 4>{x) is an “integration constant” that will be related further ahead to the initial 
entropy state. 


Evidently, the entropy integral (29) only converges asymptotically to a finite 


terminal equilibrium or “saturation” value for all fundamental observers in voids whose 
FLRW background is ACDM (or any FLRW model with A > 0). For fundamental 


observers in voids whose background has A = 0 the logarithmic decay in (|28|) leads to a 
diverging asymptotic limit I 




(In T/T)dt ~ [In r] 


—)■ cx) as r —)■ cx). 


(30) 


These results coincide qualitatively with the perturbative analysis of [10] (though an 
open FLRW background was not considered in this reference). 


II In we concluded, mistakenly, that the time integral of Sg^ converges asymptotically for hyperbolic 
LTB models with A = 0 because Sg^ —>■ 0 holds as t —)■ oo. The fact that we now show that Sgr diverges 
logarithmically for voids in an open FLRW background corrects this mistake, as these voids are based 
on such hyperbolic models. 
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Considering the radial component of dSg^ given by (14) and (19), it was shown in section 


9 of [5] that [ds^Jdxlr —)• 0 as x —t oo holds for all LTB models with A = 0 along radial 
rays: i.e. cnrves with (r, 'd, (p) constant, which are spacelike geodesics of the LTB metric 
and integral cnrves of the spacelike nnit vector n°' = orthogonal to and to the 

orbits of SO(3). 

Following |TB], it is straightforward to extend the above mentioned proof to the 
case of LTB models with A > 0 admitting a FLRW backgronnd in the asymptotic radial 
direction, as the metric fnnction a, the covariant scalars p, /C, EL and their associated 
q-scalars tend as x cxd to their backgronnd valnes a, p, H, K. (an over-bar denotes 
FLRW qnantities), while the flnctnations vanish in this limit. Considering the 
polynomial asymptotic forms in eqnations (84)-(86) and (88) of |S] together with (|^- 
(|^ we have in the asymptotic range of radial rays in arbitrary slices of constant r 

5 + i?ixV 6'^'= ^ ~-Ef-X-'' (31) 


ELq ~ 




where i/ > 1 is a constant, H = H{t) = a/a is the backgronnd Hnbble scalar and 
Hi = Hi{t), leading to 


dsgr 


d 

-K ’ 

[ d'A J 

T=const. 

La dx \ 


z/(z/ - 1) Hf 2 ^ 
X. 


H 


(32) 


Proceeding along the lines of (29), we dehne the radial entropy integral as the line 
integral of the one-form dSg 


Sgr(A)|r = 


along the radial rays: 

(dsg, ■ n) dx = 


r 

dsgr 

L 

. dx . 




(33) 


T—const 


where [dSgr/dx]T=const is given by ( [l^ and r > 0 is now an arbitrary hxed hnite 
parameter and we set the “integration constant” Sg,(r, 0) = Sgr(r)|^=o to zero (becanse 
the flnctnations vanish at the symmetry centre). Evidently, for the polynomial 

decay considered in [S] the integral Sg,.(x)|r above converges to a hnite asymptotic valne 
valne as x oo for radial rays in all FLRW backgronnds. This is a generic behavionr 
for a wide range of initial conditions and wonld only change if the latter are chosen snch 
that a and the covariant scalars decay logarithmically along radial rays (which merely 
places the restriction to avoid snch initial conditions). 


6.3. The Entropy line integral. 


The initial state of the GET entropy (29) for fnndamental observers is the GET entropy 
evalnated at Rs for all x- Hence, we can identify the integration constant in (29) with 


this initial state given by the radial integral (33) for r = 0 (or t = Hg) 

r Ids 1 

S.,U.\)=S.,(x)lo= / dx. 


(34) 


r=0 


Therefore, for any given fnndamental observer at cosmic time t = t* > ti^s (r = r* > 0) 
comoving along a worldline marked by arbitrary hxed x = A* ^ 0) fhe gravitational 
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entropy can be defined as the line integral of the GET entropy one-form dSg^ along two 
paths: a radial ray at r = 0 running from X = 0 fo X = X* the comoving worldline 
marked by x = X* running from r = 0 to r*, which is the sum of the two entropy line 
integrals (34) and ( [^ 

5i.(t*,x*) = Vi.s(x*) + v(t‘)|x*- 


(36) 


Evidently, if we are interested in the GET entropy referred to fundamental observers the 
radial entropy integral (33) only contributes to the evaluation of the total GET entropy 
(35) in the specification of a different initial entropy state (34) for each observer. As 
shown in the previous section, this radial integral generically converges, and thus it 
merely adds a finite r-dependent initial condition term to the time integral Sg,-(r)|^, 
which will not affect the qualitative asymptotic behaviour of the GET entropy for 
fundamental observers that follows from (27)-(28) and (29). In fact, by assuming 
fundamental observers in bounded comoving domains that are radially asymptotic to 
an FLRW background, we can state the following general result: 


• Sgr converges to a finite value for observers in voids with a AGDM background, 

• Sgr diverges logarithmically for observers in voids with a A = 0 background. 

which for a wide range of LTB configurations coincides qualitatively with the results of 
the perturbative analysis of HD]. 


7. Numerical analysis. 

The asymptotic results of the previous sections reveal that the growth of the GET 
entropy for fundamental observers is suppressed for generic models, but provide no 
information on the cosmic times in which this suppression sets in for interesting dust 
inhomogeneities. To look at these issues properly we need to obtain more detailed 
information on ds^JOr and its time integral Sgr(r)|;^ in their full time evolution from 
suitable initial data at ^ls (or t = 0). 


7.1. Evolution equations. 


In the case A = 0 and in some cases with A > 0 (see [23]) the form of ds^Jdr in (18) 
can be found analytically from the solutions of the Friedman equation ([^ given by the 
quadrature (21), as the latter are expressible as elementary functions that can easily be 
parametrised by the variables Ag, [71 [TH], HH] EO] EH ES El]. In the general case 
A > 0, Kqi^s 7^ 0, this task is better undertaken by means of numerical or semi-analytic 


^ Our integration of the GET entropy one-form dsgr is well posed and self-consistent regardless of 
general integrability considerations, as any one-form can be consistently integrated as a line integral 
along specific paths, even if it is not closed or exact. For the study of the frame dependent GET 
entropy on LTB models the physically meaningful paths are: radial rays at the initial time slice 
(spacelike geodesics and integral curves of n“) and the comoving worldlines (timelike geodesics and 
integral curves of u“). Discussing the general integrability of the GET entropy one-form dsgr is beyond 
the scope of this article and will be examined elsewhere. 
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solutions of ([^ (see for example m)- However, it is far more efficient for numerical work 
to compute all quantities directly (through the by the numerical solutions of 

the following system of autonomous ODE’s mQ- 

dflq 


dr 

dhq 

dr 

dr 

da 

dV 


3/ig hg. 

(36 a) 

— hq — p.q + X, 

(366) 


(36c) 

- (2hg + 3D('^))D('*), 

(36 d) 

a hq^ 

(36e) 


(36/) 

i dimensionless variables {/ig, hq, 

D(^), and the dimensionless time r are 


dehned in ( p!^ and (20). 


1.2. Initial conditions. 


The system (36a)-(36/), which is subjected to the algebraic constraints (§ , must be 


integrated numerically for initial conditions (15a)-(15&) and (16a)-(16c) specihed at 


^Ls- Since our aim is to study the GET entropy growth and its time integral for 50-100 
Mpc voids at present cosmic time tq, we need to explore which parameters of these 
initial conditions determine the entropy growth in such present day structures. For this 


purpose, we remark from looking at the asymptotic expansions (27)-(28) that the late 
time rate of GET entropy growth is basically determined by the amplitude (~ |A;i|) 
of spatial curvature initial fluctuation Dls"^ and that this growth is insensitive to the 
amplitude (~ |mi|) of density fluctuation However, this sensitivity to initial 

spatial curvature fluctuations is different when A is zero and nonzero: 

• If A > 0 a larger amplitude of Dls^ (larger |fci|) yields larger entropy growth, 

• If A = 0 larger entropy growth occurs for smaller amplitude of (smaller |fci|), 
or equivalently, for smaller deviation from spatial flatness. 

As we show further ahead, these patterns are conhrmed by the numerical examples. 
In particular, ( [2^ implies that entropy growth should be much smaller in voids in an 
open FLRW background (negative spatial curvature) than in the spatially flat Einstein 
de Sitter background, as curvature fluctuations in the latter are expected to be much 
larger than in the former. 

Since, as follows from numerical trials, different shapes of the admissible functions 
g(m) g(A:) (15a) do not produce signihcant qualitative changes in the shape of 


The evolution equations used in are equivalent to (36a)~(36d), but were constructed with the 


Aq and exact dimensionless relative functuations (called “perturbations”) 
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[^Sgr/^Tjr, we can explore effectively the effects of initial conditions to obtain entropy 
growth rates for specihc void models by varying ki for a hxed mi, and nsing the simplest 
fnnctional form for In particnlar, we select the following simple polynomial 

prohles for these fnnctions: 


VL 


(m) 


k'LS — 


+ 


mi 


1 + {x/m2y 


^LS — + 


ki 


(37) 


l + (^/fc2)3/2’ 

where the dimensionless comoving radins y is dehned in ( |T 7 | . Since the valne of 10 
an appropriate npper bonnd for the amplitnde of CDM density flnctnations at Gs when 
nsing LTB models to describe voids in strnctnre formation scenarios [IB], we have fixed 


^ is 


the maximal valne |mi| ~ 10“^, leaving the valne of |fci| as a free parameter in the initial 
conditions to be varied in order to obtain voids with ~ —0.8 density contrast of a typical 

Nnmerical trials nsing the ansatz (37) with |mi| ~ 10“^ 

10-3 


size of 50 — 100 Mpc at Tq 
reveal that the range of valnes |A;i| ~ lO-^ — 10“^ leads to snch voids. Hence, we chose 
for setting np specihc nnmerical examples the constant parameters mi, m 2 and ki, ^2 
snmmarised in Table [T] that yield the desired voids in an Einstein-de-Sitter backgronnd 
(EdS), an FLRW backgronnd with negative cnrvatnre and A = 0 (open FLRW), and a 
ACDM backgronnd. The density contrast, Hnbble scalar and mass distribntion for voids 
with these parameters are depicted by Fig. [Tj For the initial conditions snmmarised 


(a) (b) (c) 




Figure 1. The CDM void examples. The figure depicts the density contrast 
(panel (a)), Hubble scalar Ho (panel(b)) and mass distribution M = (47r/3)pqol?o 
solar mass units for the void models that follow from the initial conditions listed in 
Table 1. All curves are evaluated at present cosmic time tq as functions of area distance 
Ro = oo r. The black, red and blue curves (online version) denote the voids in ACDM, 
EdS and open FLRW backgrounds. Notice that the three void models have practically 
the same radial distribution of CDM density and mass. However, the Hubble scalar is 
different for each void, as its background value (see Table 1) has been selected to meet 
cosmic age constraints. 

in Table 1, the Big Bang time function t^,^,{r) is not simultaneous, having a background 
value Rb = Rb(oo) = 0 that corresponds to r = —2/3 and a central value tbb(O) that 
corresponds to r = —2/3 + > —2/3 (regularity conditions require t^^^ < 0), with 
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Figure 2. Initial entropy state. The figure displays the graph (solid curve) of 
the logarithm of the initial entropy state integral SgrLS defined in (38), and used to 
define the normalised GET entropy growth and its integral in (39). For the purpose 
of comparison we included the logarithmic graph (broken curve) of the initial value of 
the time derivative ds^^/dr. Notice that the latter is about three orders of magnitude 
larger than SgrLS- 



CDM prohle 

Curvature prohle 

Present background 

EdS 

aiy = 1.0 

mi = —0.005 
m2 = 0.001 

= -0.009 

k2 = 0.01 

= 1 

Hbo = 47 km/s/Mpc 

OCDM 

= 0.988250 

mi = —0.005 
m2 = 0.001 

= -0.00175 
k^ = -0.009 
k2 = 0.01 

= 0.27 

= 0 

Hbo = 59 km/s/Mpc 

ACDM 

cf?’ = 1.0 

mi = —0.005 
m2 = 0.001 

= 0 

k^ = -0.009 
k2 = 0.01 

= 0.27 

= 0.73 

Hbo = 69 km/s/Mpc 


Table 1. Parameters for the various illustrative models. These parameters 
define initial conditions at t = ^ls {r = 0) (through the ansatz ( [^ ) to obtain the deep 
voids of ^ 50 Mpc extension at present day cosmic time r = tq whose characteristic 
features are displayed in Fig. The numerical values were chosen to comply with 
|Dls^I ~ |Dls^I (or \mi\ « |fci| ~ 10“^). These are also the parameters used to define 
initial conditions in Fig. and Fig. Note that H^q (background Hubble scalar at 
t = to) is not the local Hq: it must be determined by demanding that the background 
complies with the cosmic age constraint to ~ 13.7 Gyr for values Ho"*\ supported 
by observations. 


(5tbb > 0 of the same order of magnitude as the maximal amplitude of the fluctuations: 
i.e. 5ibb ~ 0(max(|mi|, |A:i|)). Hence, the Big Bang non-simultaneity is of the order of 
to 10 “^fLS 5 or ~ 10^ — 10^ years, which is negligible in comparison with cosmic 
age [5l|20]. Since LTB models are not valid for t < fbs (or r < 0) because of the non- 
negligible presence of radiation, we only consider the evolution for t > ^ls (or r > 0), 
and thus our numerical study is not affected by dominant decaying modes very close to 
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the Big Bang time, as the latter are already subdominant at ^ls (see numerical examples 
in [3120]). 


7.3. Numerical results. 

To examine numerically the evolution of the GET entropy for arbitrary fundamental 
observers in the desired void structures (see table 1 and Fig. [^, we need to evaluate 


numerically the initial entropy state at Gg in (|34|): 

[^gLS 


■SgrLS — 


r(x) |r=0 


a 


3 

2 




+ 3D["s^],^ 


X 




.0 [V + 3D;^)],s |Vs + 3D[t^||D 


LS 


■dx, 


(38) 


for these void structures. This integral is displayed in Fig. |^for the void in the ACDM 
background, together with [c}Sg,./9r]Lg (these plots are practically the same for voids in 
the other backgrounds). 


The next step is to examine numerically the rate of growth in (18) and the entropy 


integral in (|35|) normalised by the initial state (38) 

dr. 


dS.. 


dr 




r 

'ds,; 

L 

. _ 




(39) 

J X -^grLS *grLS 

In Fig. H we display the graph of the growth rate dS^Jdr above as a function of (r, x) 
for the void structures under consideration, while Fig. depicts this growth rate and 
its time integral S'g,.(r)|^ for a typical fundamental observer (y fixed) inside these voids. 
It is clear from Fig. and Fig. that the GET entropy increases for fundamental 


ACDM 


Open FLRW 


EdS 



Figure 3. Evolution of the rate of growth of the gravitational GET entropy dS^^/dr 
as a function of (r, y) for fundamental observers inside CDM voids in a ACDM 
background, an open ELRW background and an Einstein-de-Sitter (EdS) background. 
The parameters for each case are given in Tab.[2 The red curves correspond to present 
cosmic time t = tq (or t = Iq = 13.7. Gyr). 

observers inside voids in the three backgrounds. The growth of dS^Jdr and S'gr(r)|^ is 
more intense for the void in an EdS background {dS^Jdr eventually starts decreasing for 
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(a) (b) 




Figure 4. Gravitational GET entropy for a typical comoving observer. The 

curves correspond to a fixed comoving coordinate y = 0.15 that represents an observer 
in the transition from the void centre to the background. Black, red and blue curves 
(online version) depict curves for voids in a ACDM, Einstein de Sitter and open FLRW 
backgrounds. Panel (a) displays the growth rate dSgJdr, while panel (b) depicts 
log^o ‘5'gr, both quantities defined in (39). Notice how the growth rate starts decreasing 
well before present cosmic time tq for the voids in ACDM and open FLRW background 
(it starts decreasing at r ~ 10 tq for the void with an EdS background). In panel (b) 
the logarithmic growth of is evident from the curves for the voids in the EdS and 
open FLRW backgrounds. Notice that at tq the GET entropy in the void in an open 
FLRW background is smaller than in the void in a ACDM background. However, for 
T ^ 2 to the CET entropy in the latter reaches a terminal value, while it keeps growing 
in the former. 


much larger times ~ 10 tq). This is an expected result from the asymptotic expansion 
(28), but can also be understood from the dependence of this growth rate on the Hubble 
scalar fluctuation which (from (|^) is closely connected with the radial gradients 

'H'g and 77'. Since Tiq and 77 decay at a much slower rate for the EdS background and 
(at the same time) will take large values inside the void, these gradients are necessarily 
much larger for voids in this background than in those with a ACDM and open FLRW 
backgrounds. The larger growth of dS^Jdr and for an EdS background is also 
reported in the hgures of [10] . 

It is straightforward to see in Fig. and in the panel (a) of Fig. that dS^Jdr 
in the void in a ACDM background has already begun to decrease at Tq (roughly as in 
hgure 3 of [IQ]). Remarkably, as shown by the same hgures, the negative curvature of 
the open FLRW background also leads to a suppression of the growth rate before Tq 
( this case was not examined in [ID]). As expected from the asymptotic expansions (27) 


and (28), the graphs of Fig. [^and Fig. |^con£rm that dS^Jdr —)■ 0 holds for all voids, 
thongh at very different time rates. Therefore, as shown by panel (b) of Fig. only 
for the void in a ACDM backgronnd the exponential decay is snfhciently rapid to allow 
for to reach terminal “saturation” valnes at (ronghly) r ~ 2tq for all fnndamental 
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observers. As shown in Fig. the satnration values of are position dependent inside 
the void, ranging between 10® to almost 10® times its initial state Sg,.Ls. 

For the void in an open FLRW background the logarithmic decay of dS^^/ c?r —)• 0 
leads to an asymptotic logarithmic growth of (panel (b) of Fig. |^. Hence, negative 
spatial curvature does yield a stronger suppression of the entropy rate of growth in 
comparison with the spatially flat EdS background, but it is not sufficiently strong to 
yield a finite terminal value for The time prohle of the curves of panel (b) of Fig. 
for the ACDM and EdS backgrounds are qualitatively analogous to the first and fourth 
curves (top and bottom) of figure 2 of [10] that correspond to these backgrounds. 


8. Comparison with the perturbative treatment. 


It is important to compare the results obtained in previous sections with those of the 
perturbative analysis in [101 . If we bear in mind the correspondence between fluctuations 
of the q-scalars and spherical dust perturbations of the linear theory (see detail in 
[T7]h the common theoretical ground (related to the Weyl tensor) between the non- 
perturbative and perturbative results can be appreciated by a direct comparison between 
the conditions for entropy growth in both treatments. Since the Weyl tensor and its 
squared contraction for LTB models take the form [5] 




/l[“ — ?tU\rV}°' 








it is then evident that equation (10) (and similar equations in [SIE]) and equation (37) 


of [To] are equivalent at first order 


k(|D<qan) 



(41) 


where a and "H are the FLRW scale factor and Hubble scalar, p = J dt/a is the conformal 
time and we used the fact that at linear first order we must have (see ini) : a ~ a, F ~ 1, 
as well as {CahcdC°-^‘^‘^)v ~ CabcdC°'^^'^ cc where the overline in ( )x) denotes an 

ensemble average (see [T0]L 

The connection between the perturbative and non-perturbative approach to the 
GET entropy growth can also be appreciated from the relation between the exact LTB 
evolution equations and evolution equations for linear perturbations. To highlight this 


relation we recast the system (36a)-(36d) into a dimensional form that is suitable for 


comparison with evolution equations for dust perturbations in a synchronous comoving 
gauge (see comprehensive discussion in di): 

t = (42a) 

— — A-it 

= -nl - -p.. + yA, (426) 

=-3(1+ (42c) 

f cnn An I I t\{H) o(nj nj \2 dn _ 
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where p^^it), Tilsit) are the limits of p, Ti (and also of Pq, T-iq) as r —)■ cx) for arbitrary- 
fixed f, and thus correspond to the density and Hubble scalar of the asymptotic FLRW 
background 


(in fact, (42a)~(426) determine this background), while 5^Jf\ are the 
asymptotic non-local fluctuations defined as 
^(p) ^ P — Pas 




Pa 


(43) 


whose direct connection to the gauge invariant density and Hubble scalar perturbations 
of linear perturbation theory is discussed extensively in im (in particular, also 
corresponds to the density contrast). 

A linear regime associated with (42a)-(42d) and (43) can be defined by the 
conditions <C 1 and |DiI^V^as| ^ 1, which imply (as long as they remain valid) 
that the quadratic terms and are negligible, and thus H ^ Hq ^ 

hold. Hence, the dynamics of the dust fluctuation (in this regime) is determined by the 
linearised forms of (42c)-(42d): 

5^p'> = 


3D(«) i)^ =-2 KsD(^)-(44) 

as * as as ^ r as ' \ / 

which, as shown in ini. are fully equivalent to the evolution equations of standard dust 
cosmological perturbations in the isochronous comoving gauge. In fact, if combined (44) 
lead to the well known second order equation — Anp^S^P^ = 0 for linear 

dust perturbations. 

Since void structures defined by the initial conditions at ^ls listed in Table 1 evolve 
into a fully non-linear regime at tq (see Fig. [^, it is straightforward to show that these 
structures necessarily comply with a linear regime for restricted evolution times that 
are close to the initial time Gs (or r ~ 0). Hence, we can state unequivocally that the 
entropy growth rate, dS^Jdr, and the entropy integral, 5gr(T)|^, reach much smaller 
values in the linear regime than in the non-linear regime. This becomes evident by 
comparing the early times (r ~ 0) linear regime values of these quantities with their 
non-linear regime latter times (r > ro/2) in their plots of Fig. |^and Fig. 

However, a more accurate quantitative comparison between the linear and non¬ 
linear regime can be achieved by looking at the evolution of the GET entropy for void 


structures that remain in the linear regime up to tq. From the discussion of section 7.2 


it is reasonable to guess that the right initial conditions for such voids involve choosing 

(k) 

values of the amplitude \ki\ of the initial spatial curvature fluctuation Dls^ that are 
smaller than those of Table 1, since this amplitude is the determinant parameter for the 
rate of growth of the GET entropy. 

Gonsidering now only voids in a AGDM background, we build up initial conditions 


IU'ls I 


< 




as described in the previous paragraph: fluctuations complying with 
0(10“^), which can be achieved by the ansatz (37) with \ki\ < \mi\ = 0.001. For such 
initial data we solve numerically the linear evolution equations (42a), (42d) and (44). 
The density contrast and Hubble scalar of the resulting void configurations are displayed 
by Fig. i where the curves marked by B, G and D correspond to the following decreasing 
values: ki = —9 x 10“^, —9 x 10“®, —9 x 10“®, with the parameters mi, m 2 , ^2 kept the 
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same as the void in the ACDM background in Table 1. For the purpose of comparison, 
we include the curves marked by A that correspond to the non-linear void in a ACDM 
background used in Figs. Bi and 1^ As shown by the curves B, C and D of Fig. the 
voids remain linear (he. very “shallow”) up to present cosmic time tq. 

In order to compare the GET entropy growth of the voids in the linear and non¬ 
linear regimes, we plotted in Fig. [^the saturation terminal value of (panel (a)) and 
the saturation per comoving volume (panel (b)) for fundamental observers in the voids 
of Fig. Since reaches its saturation at r ~ 2ro, the saturation values are obtained 
by evaluating at r = Stq, whereas the ratio of saturation to comoving volume is 
obtained by plotting S'g,(3ro)/i?^(3ro), where i? = ar is the area distance. As expected, 
the GET entropy growth in linear voids is much smaller than in non-linear ones, as the 
magnitudes of S'gr(3ro) and S'gr(3ro)/i?^(3ro) for the curves B, C and D in Fig. [^decrease 
(with respect to the curve A) in direct proportionality to the decrease of the amplitude 
parameter \ki\. 

In particular, it is important to emphasise that the numerical results of panel (b) 
can be directly compared to the numerical results of [lO] , who only examined saturation 
GET entropy values per comoving volume. Notice how the curves B, G and D in panel 
(b) of Fig. 1^ for the linear voids rapidly decrease to values of ~ 10“^, which is the 
numerical value displayed in the fourth curve (top to bottom) of hgure 2 of [lO] for the 
GET saturation per comoving volume in linear perturbations in a AGDM background. 
Notice that also in the non-linear void (curve A) the saturation values per volume 
decrease to the perturbative values ~ 10“^ for increasing comoving volume, as the non¬ 
linear void becomes linear at large scales where its parameters approach the background 
parameters. 

As a test of consistency, we verihed that exactly the same curves of Fig. and Fig. 
are obtained (for same initial conditions) from the full exact evolution equations (either 
(36a)-(36d) or (42a)-(42d)). We also verihed (but are not displaying the graphs) that 
Agr diverges logarithmically for linear voids in an open FLRW and Einstein de Sitter 
backgrounds. Hence, our numerical results fully agree with those of [TO] . 


9. Conclusion. 

We have examined the exact non-perturbative effect of negative curvature and A 
term on the evolution of the GET gravitational entropy, Sg^, inside of expanding sub¬ 
horizon 50-100 Mpc voids that form in three prototypical standard FLRW cosmological 
backgrounds: AGDM concordance-type model, open FLRW (negative spatial curvature 
but A = 0) and Einstein de Sitter (EdS) (see Table 1 and Fig. [^. We used analytic 
expressions to examine the asymptotic time behaviour of Sg,. and its line integral Sg,.(f)|r, 
and a numerical study to determine their full evolution from linear initial conditions 
specihed at the last scattering time t = Gs- 
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(a) (b) 




Figure 5. Voids in the linear regime. The figure depicts the same quantities of 
panels (a) and (b) of Fig. (density contrast and Hubble scalar), but for voids in a 
linear regime and in a ACDM background (the curve marked by A corresponds to the 
non-linear ACDM void used in previous figures). The curves marked by B, C and D 
correspond to the linear voids whose parameters are defined in the text. Notice that 
the voids associated with the curves B, C and D have the negligible density contrast 
and radial variation of "Ho expected from the linear regime. 


(a) (b) 




Figure 6. Saturation values. Panel (a) depicts the logarithm of the terminal GET 
entropy saturation values at t = Stq) for fundamental observers in the voids of 
Fig. El Panel (b) displays the same saturation values weighed by their comoving volume 
defined by where R = ar is the area distance at Stq. Notice that the saturation 
values decrease several orders of magnitude for the linear regime voids (curves B, C 
and D). On the other hand, the saturation value per volume reaches for the linear 
voids (and the non-linear one at large scales) the same order of magnitude numerical 
values ~ 10“"^ reported by in Fig 2 of the perturbative study of m- 
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We have found that for all such voids the rate of entropy production is initially 
growing, but this growth has started to be suppressed before our cosmic time t = to in 
voids in the ACDM and open FLRW backgrounds, while for voids in an EdS background 
this suppression only begins at very late times t ~ 10 to- While for all voids we have 
Sgr —>■ 0 asymptotically, this decay occurs at a rapid exponential rate for voids in a 
ACDM background and it occurs at a slow logarithmic pace for voids whose FLRW 
background has A = 0 (open FLRW and EdS). As a consequence, only for voids in 
a ACDM background this asymptotic decay is sufficiently rapid to allow Sgr(t)|r to 
converge asymptotically for all fundamental observers (at about f ~ 2fo) into a terminal 
equilibrium (or “saturation”) value between 10® and 10® times its initial state defined 
by the radial entropy integral ( |^ at f = Rs fixed (see Fig. |^. For voids in the open 
FLRW and EdS backgrounds Sgr(t)|r diverges logarithmically (though at a much slower 
pace for the former). These results are all depicted graphically by Fig. |^and Fig. 

Since our results should agree in the appropriate regime with those of [10] based 
on linear dust perturbations, we provided a direct comparison between the exact non¬ 
linear evolution of the CET entropy for voids (in an ACDM background) that are fully 
non-linear at R and voids (whose initial conditions are defined also at Rg) that remain 
in a linear regime up to to (see Fig. [^. The results of this comparison are as expected: 
as shown in Fig. [^the saturation terminal values of the CET entropy decrease for the 
linear voids in the same order of magnitude proportion as the amplitude of the initial 
spatial curvature fluctuation. However, we also obtained a good quantitative agreement 
with the results of [10] : the saturation values of the CET entropy per comoving volume 
(the curves B, C and D in panel (b) of Fig. rapidly decay to values of ~ 10““^ that 
fit nicely with the same quantity plotted in curves in Fig 2 of ra (see specially the 
fourth curve from top to bottom that corresponds to the ACDM background). In fact, 
this agreement also holds for the non-linear void (curve A), which shows that linear 
perturbative conditions should also hold in the large scales for non-linear structures 
where the latter converge to the FLRW background. 

It is worth mentioning from a theoretical perspective that the asymptotic value of 
the CET gravitational temperature for voids in a ACDM background is proportional to 
the cosmological constant (this temperature is zero for voids in A = 0 backgrounds). 
This fact indicates the existence of a potential link between the value of the cosmological 
constant and thermodynamical considerations, supporting the notion that gravity and 
thermodynamics might be intertwined at a fundamental level; see e.g. ESI for attempts 
at such an emergent gravity scenario. 

Another important theoretical issue is the fundamental nature of the gravitational 
entropy and its connection to other entropy concepts in the context of gravitational 
systems, such as the holographic Hawking-Bekenstein black hole entropy [2Z1I2H] and 
its various extensions to fluid sources (in FLRW spacetimes [29] and in LTB models 
[a El EQ]). Since this issue is still an open problem whose proper treatment would 
require a separate article to be submitted in future work, we introduced in section [^ 
a fiducial fundamental constant (or scale parameter) with entropy units, cr, in order 
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to examine the evolution of the GET entropy through dimensionless expressions. As a 
consequence, the numerical values we have obtained for the GET entropy need to be 
understood as entropy values proportional to this yet unknown fundamental entropy. 
See [31] for a recent way of relating the Hawking-Bekenstein entropy of a black hole to 
the variation of the GET entropy outside the collapsing dust at the origin of the black 
hole. 

Evidently, we have only examined very idealized spherical expanding cosmic voids, 
and thus further research is needed to probe the GET proposal (and the proposal of 
[1]) on more general spacetimes, such as Szekeres models [22|, and on the process of 
structure formation and gravitational collapse. In particular, we aim at studying the 
growth of these gravitational entropies in the context of the formation of virialized 
stationary structures, which may provide a connection with theoretical work done on 
n-body numerical simulations na and Newtonian self-gravitational systems [32], as 
well as research on various proposals on non-extensive entropy definitions [33]. This 
research is currently under way and will be the subject of future works. 
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